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We study the 1 + 1-dimensional random directed polymer problem, i.e., an elastic string (f>(x) 
subject to a Gaussian random potential V((j), x) and confined within a plane. We mainly concentrate 
on the short-scale and finite-temperature behavior of this problem described by a short- but finite- 
ranged disorder correlator U((j>) and introduce two types of approximations amenable to exact 
solutions. Expanding the disorder potential V(<j>, x) « Vo(x) + f(x)<j>(x) at short distances, we 
study the random force (or Larkin) problem with Vo(x) = as well as the shifted random force 
problem including the random offset Vo(x); as such, these models remain well defined at all scales. 
Alternatively, we analyze the harmonic approximation to the correlator U(<j>) in a consistent manner. 
Using direct averaging as well as the replica technique, we derive the distribution functions Vl, v (F) 
and Vl(F) of free energies F of a polymer of length L for both fixed (</>(£) = y) and free boundary 
conditions on the displacement field (j>(x) and determine the mean displacement correlators on 
the distance L. The inconsistencies encountered in the analysis of the harmonic approximation 
to the correlator are traced back to its non-spectral correlator; we discuss how to implement this 
approximation in a proper way and present a general criterion for physically admissible disorder 
correlators U((f>). 

PACS numbers: 05.20.-y 75.10.Nr 74.25.Wx 61.41,+c 
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I. INTRODUCTION 

Directed polymers subject to a quenched random po- 
tential have been the subject of intense investigations 
during the past two decades^. Diverse physical sys- 
tems such as domain walls in magnetic films 2 -, vortices in 
superconductors^, wetting fronts on planar systems^, or 
Burgers turbulence^ can be mapped to this model, which 
exhibits numerous non-trivial features deriving from the 
interplay between elasticity and disorder. The best un- 
derstanding, so far, has been reached for the (1 + 1)- 
dimensional case, i.e., a string confined to a plane, and it 
is this geometry we study in the present paper. Specif- 
ically, we analyze the situation illustrated in Fig. [T] an 
elastic string (with elasticity c) of finite length L within 
an interval [0,L] directed along the x-axis. The disor- 
der potential V(4>, x) drives a finite displacement field 
4>(x), which is counteracted by the elastic energy density 
c(d x 4>) 2 /2. The problem is conveniently defined through 
its Hamiltonian 



H[<j>(x);V} = I dx{^[d x <fi(x)] 2 + V[cj>(x),x]}; (1) 

the disorder potential V(cf), x) is Gaussian distributed 
with a zero mean V(cj>, x) = and a correlator 



V(4>, x)V(cf>',x') = 5{x - x')U{4> - (/>'), 



(2) 



with U((j>) the correlation function. In the present work, 
we are mainly interested in short-range correlated disor- 
der potentials, which we characterize by their extension 
£ and the strength Uq = U(0); these parameters then 
combine in a curvature u = —U"(0) « f/(0)/£ 2 . 




FIG. 1: Thermally averaged trajectory (<j)(x))th of a random 
directed polymer in a fixed disorder potential V(4>, x) start- 
ing in (x,(f>) — (0,0) and ending in (L,y). The free energy 
associated with such a configuration is denoted by F. The 
random choice of the underlying disorder potential V(cj>, x) 
defines a random process; the free energy then turns into a 
random variable, whose distribution function VL, y {F) we seek 
to calculate. 



Quantities of central interest a re th e scaling behavior 
of the mean-squared displacement {<p 2 ){L) of the polymer 
with length L and the distribution function P(F) of the 
polymer's free energy F. For free boundary conditions 
at x — 0, L, the former is defined through the expression 
W)(L) = ([0(£)- </>(0)] 2 ), where (...) and {777) denote 
thermal (temperature T) and disorder (random potential 
V) averages, respectively. Its dependence on the distance 
L exhibits a scaling form (<fi 2 )(L) oc L 2 ^, with the so- 
called wandering exponent £ to be determined. 6 

The polymer's free energy F is defined via its partition 
function 

fi>(L)=y 
Z[L,y;V}= V\<j>(x)] exp^tf^*); V}), (3) 

J0(O)=O 

where /3 = 1/T denotes the inverse temperature (we set 
the Boltzmann constant to unity), from which the free 



energy 



F[L,y;V] = ~-T\n(Z{L,y;V}) 



(4) 



follows immediately. The free energy F in Eq. Q is 
denned for a specific realization of the random potential 
V and thus defines a random variable; given the above 
(Gaussian) distributed disorder potential, the task then 
is to determine the distribution function VL,y{F). In Eq. 
([3]), we have considered a string starting at (x, <p) — (0, 0) 
and ending in a fixed position (x, <f) — (L, y) a distance 
L away, but other cases, e.g., a free boundary condition 
at x — L, see below, may be studied. 

Two types of analytic solutions are known for the (I + 
I)-dimensional random polymer problem with a short- 
range, in fact, (S-correlated, potential U{<j>) — u8{4>): 
mapping the replicated problem to interacting quantum 
bosons^ and using the Bethe-Ansatz technique, one can 
find the spectrum and eigenfunctions of the interacting 
quantum many-body problem, from which the distribu- 
tion function Vl{F) for the free energy F of a polymer 
of length L and fixed endpoint at <j)(L) = can be ob- 
tained; we call this the 'longitudinal problem'. Restrict- 
ing the solution to the ground state wave function only 
permits the determination of the far- left tail 8 . First in- 
dications that the full distribution function should be of 
the Tracy- Widom form derived from the work of Prahofer 
and Spohn 9 - on polynuclear growth, a model in the uni- 
versality class of the Kardar-Parisi-Zhang (KPZ) model, 
to which the random directed polymer problem belongs 
as well. Recently, both tails of the free-energy distri- 
bution function Vl{F) have been found using instanton 
techniques 10 , with results consistent with those in Ref. 
9. Finally, the full distribution function (in the Tracy- 
Widom form) has been obtained independently by sev- 
eral groups, by Dotsenkc 11 and by Calabrese et al~2. us- 
ing the replica technique including the full spectrum, by 
Sasamoto and Spohn 13 deriving the exact KPZ solution 
from the corner growth model, and by Amir et ali*. 

An alternative (exact) result has been obtained using 
a mapping to the Burgers equation via the Colc-Hopf 
transformation 15 ; making use of an invariant distribu- 
tion, the distribution function V y (F') for the free-energy 
difference F' between two configurations with endpoints 
separated by 2y has been foun d 15 ' 16 (the so-called 'trans- 
verse' problem). Both approaches have been helpful in 
finding the wandering exponent^ 5 - £ = 2/3 of transverse 
fluctuations of the polymer; it is generally believed that 
this value is universal for short-range correlated disorder 
potentials, i.e., for rapidly decaying correlator functions 
U(\<f>\ -+ oo) -► 0. 

Recently, some of usl£ have studied the joint dis- 
tribution function VL,y(F,F') for a polymer of length 
L which involves two configurations of the string end- 
ing in points separated by 2y. An interesting result is 
found for the ^-correlated potential, where Vl, v {F, F') — 
Pl, v (F) Py(F') separates for large L and large negative 
values of the mean free energy F, with the factor Pl^ v (F) 
reproducing Zhang's tail 8 - for Vl{F) and the transverse 



part Py(F') coinciding with the stationary distribution 
function V V {F') of the Burgers problem^ 5 -, thus placing 
previously known results into a common context. 

The above solutions focus on the (^-correlated disor- 
der potential which exhibits a singular zero-temperature 
limit; the problem with a finite-width correlator U((j>) of 
the random potential has remained unsolved so far. In 
the present paper, we study another problem, originally 
proposed by Larkin 19 , that is amenable to a complete 
exact solution. This problem deals with the polymer's 
behavior at short distances and also exhibits a regular 
low-temperature limit. The basic idea is to linearize the 
problem, either by an expansion of the original random 
potentially V(<f>, x) for small values of <fi, 



V(cf>,x)ttV (x) + f(x)4>(x), 



(5) 



or by an expansion of the correlator U{4>), see below. 
Studying the destruction of long-range order due to the 
presence of quenched disorder, i.e., the behavior of the 
displacement field <p(x) at large distances x, the random 
shift Vq (x) can be dropped and one arrives at the Larkin 
or random-force model described by the Hamiltonian ([!]) 
with 

Vfax) = f(x)<f>(x), (6) 

where f(x) denotes a (Gaussian) random force field with 
zero mean f(x) = and a correlator 



f(x)f(x') = u8{x-x). 



(7) 



Its free-energy distribution function V Vi l(F) has been 
calculated by Gorokhov and Blatter— for fixed bound- 
ary conditions; here, we extend this analysis to describe 
the case of free boundary conditions. Furthermore, in 
our analysis of the Larkin model as a short-scale approx- 
imation for the random potential problem, we are not 
allowed to ignore the random shift Vo(x); in our study, 
we assume the latter to be Gaussian correlated, 



Vo(x)V (x>) = U 8(x-x'), 



(8) 



and uncorrelated with the force, Vo(x)f(x') = 0; we call 
this approximation the shifted random force model. Both 
models not only serve as approximations to the random 
polymer problem on short scales but also describe phys- 
ical problems where the underlying randomness is de- 
scribed by a (shifted) random force field on all length 
scales. Therefore, below we will quote results for these 
models for arbitrary lengths L. 

Another approximation discussed in the literature!^ 
deals with the expansion of the correlator U{<j)—<j) 1 ) rather 
than the potential, 

U(<j> - <f>') ~ U p {cj> - 0') = J7 - \u{4> - 0') 2 ; (9) 

also, models have been studied with ad hoc power-law 
expressions for the correlator—. Here we point out, that 



expanding the correlator U(4>) of the random potential or 
choosing arbitrary forms for the correlator is a problem- 
atic step, as this action may generate a non-spectral cor- 
relator introducing an ill-defined Gaussian measure and 
thus lead to an unphysical model at the very start, cf. 
Sec. IIVI below for a detailed discussion. As an imme- 
diate consequence, the mapping to the quantum boson 
problem requiring integration over the random potential 
V(4>, x) fails. On the other hand, performing the inte- 
gration over the random potential V{4>, x) with a well 
defined correlator U(4>) first and expanding the resulting 
interacting —f3 2 U(cf>) between bosons only thereafter is a 
perfectly admissible program producing identical results; 
in this case, however, we know that the (perfectly well de- 
fined) quantum boson problem does not describe the ran- 
dom polymer problem on scales where the approximate 
quadratic correlator deviates strongly from the original 
correlator. 

In the following, we will discuss the harmonic- 
correlator approximation of the original random poly- 
mer problem with the understanding, that the har- 
monic approximation is done after the mapping to 
bosons. This approximation then relates to the 
above shifted random force model with a force cor- 
relator c^(0, x) dp V (<!>', z')U,*'=o = [-J7"(0)] S(x - 
x') = u5(x — x') and the shift described by 



V{4>, x)V((j)', x'), j /=0 = Uq 5(x — x 1 ) with the additional 
advantage to preserve the translation invariance of the 
problem (note, that the shifted random force model in- 
volves only the constant and mixed (u<fi(j>') terms in the 
correlator, with the quadratic terms ex (<j> 2 + (/)') ab- 
sent). Furthermore, in Sec. IIV( we will present a cri- 
terion assuring the consistent definition of a correlator: 
such correlators have to be spectral. Also, we emphasize 
the difference in terminology introduced above: while the 
random force and shifted random force approximations to 
the random potential V(<fi, x) define proper models of dis- 
ordered elastic systems, this is not the case when expand- 
ing the potential correlator £/(</>); this is why we refrain 
from considering a model with a harmonically correlated 
random potential but prefer to talk about a harmonic 
approximation to the correlator. 

The (shifted) random force models and the harmonic- 
correlator approximation produce similar results for var- 
ious quantities defined at short scales, such as the mean 
free energy and the displacement correlator. While all 
displacement correlators are identical for the two models 
and the approximation, the random force model supplies 
us with the distribution function for the relaxational free 
energy (i.e., the free energy of the distorted string re- 
duced by the energy of the straight string), whereas the 
shifted random force model and the harmonic-correlator 
approximation provide the distribution function for the 
total free energy. Comparing the latter two, we note 
that the results for the harmonic-correlator approxima- 
tion provide a more consistent description of the origi- 
nal problem Eq. (J2J) at short distances, as the correlator 
remains translation invariant, while some terms in the 



expansion of the random potential are dropped for the 
shifted random force problem. On the other hand, the 
applicability of the harmonic-correlator approximation is 
limited to short scales, as the free-energy distribution 
function Vl{F) suffers from a negative second moment 
whenever large displacements show up, e.g., at large dis- 
tances L or at high temperatures T; this is to be ex- 
pected as the harmonic correlator deviates strongly from 
the original correlator U{<p) (and eventually turns nega- 
tive) at large arguments. 

For a system with quenched disorder it is usually ex- 
tremely difficult to find averaged physical quantities, e.g., 
the mean free energy F — —Tin Z . Replica theory, re- 
quiring calculation of the disorder-averaged n-th power of 
the partition function Z n then comes in as a helpful tech- 
nique. Usually, it is the limit lim n _>o [{Z n — l)/ n ] = m ^i 
to be calculated after analytic continuation of n, which 
is of fundamental interest. It turns out (see below) that 
the same quantity Z n and its analytic continuation is 
relevant in the calculation of the free-energy distribu- 
tion function V(F), since the latter is nothing but the 
inverse Laplace transform of the former—, hence replica 
theory seems the technique of choice for the solution of 
the present problem as well. However, the shifted ran- 
dom force model defines a quadratic problem that can 
be analyzed in a straightforward manner, i.e., the parti- 
tion function Z[f{x)\ (involving an integration over the 
field) can be found for any configuration f(x) of the ran- 
dom force and the disorder average of its n-th power can 
be done in the end. This is opposite to the replica ap- 
proach where the integrations are interchanged, with the 
first integration over the disorder of the replicated sys- 
tem generating an interacting imaginary-time quantum 
boson problem, which then is solved in a second step 
(corresponding to the integration over the field). Below, 
we will discuss both procedures for the Larkin model and 
find that they provide similar challenges. 

The disorder (u and Uq) and elastic (c) parameters of 
the above random polymer problems define convenient 
and physically relevant length and energy scales: The 
ratio of Uq and u defines the transverse length scale £ 
where the shifted random force model approximates well 
the random polymer problem, 



*-%r 



(10) 



Comparing the elastic energy E c = c£ 2 / 'L = cIIq/uL 
with the disorder energy Ef = \/UqL accumulated over 
a distance L, one obtains the corresponding longitudinal 
scale L Cl 



L r = 



2 C/ \i/3 /c 2 £ 2 \V 



(11) 



Finally, the energy scale associated with these length 
scales is 



U c 



(culy/^ce 

L c 



\ u 



(12) 



Note that the longitudinal (L c ) and transverse (£) scales 
define the limits of validity where our expansions describe 
the original random polymer problem. The parameters 
are not fully appropriate to describe the results of the 
Larkin model, as the latter is characterized by one disor- 
der parameter (u) only — to allow for proper comparison, 
below, we will nevertheless express all physical results 
through £, L c , and U c . For the Larkin model, these pa- 
rameters will combine to expressions containing only u 
and c. 

Besides providing new results for the Larkin model and 
a discussion of its use as an approximation to the ran- 
dom potential problem at short scales, the present study 
also has its merits from a methodological point of view, 
since this is the only case where the entire analysis (direct 
and via replica) could be carried through in a complete 
and consistent manner. Below, we introduce the formal- 
ism (Sec. |ll]) and then apply it to the (shifted) random 
force models (Sec. IIII[) . We then analyze the harmonic- 
correlator approximation (section IIV|) . analyze its failure 
due to its non-spectral property, and state the spectral 
condition to be satisfied by a properly defined random- 
potential correlator U((f>); furthermore, we briefly present 
the results for the displacement correlators which are 
idential in all three cases. Conclusions are presented in 
Sec. El 



II. METHODOLOGY 

Evaluating the partition function Eq. ([3]) and the ex- 
pression Eq. (|3]) for the free-energy F for a given random 
potential V(<f>, x) defines a sample-dependent random 
quantity. Its free-energy distribution function VL,y{F) 
can be derived from the n-th powers of the partition func- 
tion 



Z n (L, y) = Z[L, y; V) n = exp {-nj3F[L, y; V]); (13) 

these are equal to the (bilateral) Laplace transform of 
the free-energy distribution function VL,y(F) at integer 
multiples of /?,— 



Z n (L,y) 



dFV L , v (F) exp(-n/3F). 



(14) 



Hence, the inversion of this expression through the (in- 
verse) Laplace transformation provides us with the free- 
energy distribution function VL,y{F). This requires us 
to analytically continue the expression Z n (L,y) for the 
moments from integer values of n to the complex n-plane, 
Z n (L, y) — > Z(ji; L, y) with j3n — > n. The free-energy dis- 
tribution function "pL,y(F) then is given by the inverse 
Laplace transformation 



"^ = h 



-\-ioo 



dr)Z( m L,y)exp(r)F), (15) 



where the integration goes over the imaginary n-axis with 
Re(n) chosen in such a way as to place all singularities 



in Z(r/;L,y) to its right. Furthermore, taking the fc-fold 
derivatives of Z(n;L,y) with respect to ?y provides us 
with all the moments 



(F")(L,y) = (-l) 



k d k Z(r,;L,y) 



drj k 



»7=0 



(16) 



The calculation of the moments Z n (L,y) involves in- 
tegrations over the displacement field (f>(x) and over the 
distribution function P[V(tfi, x)] of the disorder potential, 



Z n (L,y) = I V[V{<t>,x)]P[V{<f>,x)\ (17) 

n » n 

J] / £[<A a (z)]e X p(-^ff[^(aO])- 



0=1 



For the random force or Larkin model, these integrations 
can be done straightforwardly in the sequence above us- 
ing the distribution function for the random force 



P[/(x)] oc exp(- fdxf(x)/2u 



(18) 



this program will be carried through in section IIII Al be- 
low. Fixed and free boundary conditions are conveniently 
imposed by the requirements 0(0) = 0, <j>(L) = y and 
0(0) =0, [d x <f>]{L)=Q. 

On the other hand, for the general situation with a 
random potential V((j), x), the integration in Eq. (|17p 
over the displacement fields cf) a (x) cannot be done. In- 
terchanging the integrations over V(4>, x) and <f> a (x) takes 
us directly to the replica technique: performing first the 
integration over the disorder potential V(<fi, x), the re- 
maining integrations over the fields 4>a(x) have to be done 
with the replica Hamiltonian H n [{4> a }], 



*({yaY,L) 



H n [{Mx)}} 



n / v\m*)\ 

x exp(-l3H n [{(t) a (x)}} 



f L { c ™ 

/ dx \ nEt^w]' 



(19) 



(20) 



J2 U[</, a (x)-if> b (x)] 



a,b=l 

Here, we have allowed the individual replicas of the elas- 
tic string to end in different locations y a . The expression 
Eq. (|19p is identical with the imaginary time (x) propaga- 
tor ty({y a }', x) of a many body problem in a path integral 
setting. Collapsing the end-points y a = y, this propaga- 
tor coincides with the n-th moment (|17|) of the partition 
function, 

V({y a =y};x = L) = Z r (n; L, y) = Z n (L, y), (21) 

where the last equation holds, provided that the inte- 
gration over the disorder potential V can be exchanged 
with the integration over the field <\> a . Note that it is this 



mapping from the polymer statistical mechanics prob- 
lem to the quantum boson problem which fails when the 
correlator is non-spectral, e.g., for the (naive version of) 
harmonic-correlator approximation. 

The equivalence to a quantum many body problem be- 
comes more obvious when going from the path-integral 
Eq. (|19l) to an operator formalism; the evaluation of the 
path-integral Eq. (fH)|) then is equivalent to the solution 
of the imaginary-time Schrodinger equation 



-d x *({y a };x) = H*({y a };x 

with the Hamiltonian 

32 



(22) 



H 



a— 1 



,6=1 



The Hamiltonian (|23|) describes n particles of mass j3c in- 
teracting via the attractive two-body potential ~f3 2 U(y); 
the propagation in (|22[) starts in the origin at time x = 0, 

(24) 



*({ 2/a };0) = n™ =1 5(y a ), 



and ends at different coordinates {y a } after propagation 
during the time x. To keep up the formal distinction be- 
tween the two quantities, we denote the direct physical 
definition of the moments by Z n (L,y) (first integration 
over the field <f>, raising the result to the power n, and 
averaging over disorder V) and denote the replica expres- 
sion (n-fold replication followed by averaging over disor- 
der V and integration over the fields <j> a done in the end) 
by Z r (n;L,y). 

Finally, we note that in the replica technique, free 
boundary conditions at the endpoint are more conve- 
niently implemented through an integration over y; the 
partition function for the polymer with free boundary 
conditions assumes the form 



+oo 



Z[L;V}=N dyZ[L,y;V] 



(25) 



with N a suitable normalization constant. Taking the n- 
th power and averaging over the disorder potential V pro- 
vides us with the moments Z n (L). Following the replica 
procedure, after replication and integration over V, one 
arrives at the free-boundary replicated partition function 
through integration over the set {y a } of n different end- 
points, 



Z r (n;L) 



\{N 



+ 00 



dy a 



*({y a };x = L). (26) 



In the next sections, we are going to apply the above gen- 
eral schemes for the calculation of the free-energy distri- 
bution functions VL,y{F) and Pl(F) for fixed and free 
boundary conditions, respectively, of the (shifted) ran- 
dom force model Eqs. (H])-© and of the random directed 
polymer model Eq. (JlJ with the parabolic approximation 
for the correlation function, Eq. ([9]) , done after averaging 
over the disorder. Before doing so, we briefly discuss the 
results for the free string which determines our normal- 
ization N . 



A. Free string 

The path integrals Eqs. ©, ([T7]). and flTUJ} over tra- 
jectories <fi(x) involve an arbitrary measure of integra- 
tion. Here, we choose a particular normalization such 
that the partition function Zo(L,y) (or wave function 
^o(y, x — L)) of the free polymer problem with fixed 
boundary conditions 0(0) = and <p(L) = y assumes the 
form 



z (L, y ) = tf (i/;£) 



(27) 



y r p c r L 2 

T>[c/)(x)] exp -— / dx[d x <j>{x)] 
o L * Jo 



= exp 



2L 



(--» 



the corresponding free energy then is given by 

C n 



F (L,y) 



21 



V 



(28) 



For the partition function of the free polymer with free 
boundary conditions we choose the normalization N — 
(2ttL/(3c) 1/2 and obtain 



Z (L) 



I2ttL 



13c J_ 



+ 00 



dyZ (L,y) = 1 



(29) 



and the free energy Fq(L) = 0. These results will be 
helpful in the interpretation of the free-energy distribu- 
tion functions for the random force model calculated be- 
low. With this normalization, all our free energies F 
are measured with respect to the free thermal energy 



F fe 



T In \j2itLT/c of the free string due to its entropy. 



III. RANDOM FORCE MODEL 

We select the simplest case, the Larkin model, for the 
discussion of the two methodological approaches involv- 
ing either direct integration over the field <j> and subse- 
quent disorder average over V or the route following the 
replica approach. While the first route is preferably done 
in Fourier space, the replica calculation will be formu- 
lated in real space. Also, note that the analysis for the 
Larkin- or random force model provides the distribution 
function for the relaxational free-energy F — Eq rather 
then the (total) free energy F of the polymer, 



Z[L,y:V] = exp"^^*^ 
r 4>{L)=y 

■D[d>(x)] ex V (-0H[(f>;Vo + f<j>}) 
/<Ko)=o 

r4>(L)=y 



(30) 



-PE 



X>foKaO]exp(-/3iJtoM 



0(O)=O 



with Eq = J dxVo (x) the disorder energy of a straight 
string. This latter remark is relevant in the comparison 
of the random force and the harmonic models. 



A further speciality of the Larkin model is the separa- 
tion between the thermal and the quenched disorder 22 . 
Indeed, splitting the displacement field 4>{x) into the 
Hamilonian's minimizer (jj q (x), 



cdl^ix) = f(x), 



(31) 



and fluctuations 54>{x) around it, 4>{x) — 4> q {x) + S(j)(x), 
we can decompose the Hamilonian into the free part 
Hq[5(/)(x)] and the energy of the minimizer H[<fi q (x)], 



H [8<P(x)} = / dx -[d x 8<j>ix)Y (32) 

Jo z 

H[Mx)] = J dxfyd^x)] 2 + fix)^)}. 



In addition, we can account for the boundary condition 
4>(L) = y through a simple shift 4>{x) — > xy/L + <j>{x), 
which adds the terms 



H v=2i+lJ dxxf(x) 



(33) 



to the Hamiltonian. The partition sum then naturally 
separates into thermal and quenched-disorder averaged 
factors , 22 



Z(L, y; /) = Z (L, y) exp{-/3(iJ[0 q (x)] + H y ) }, (34) 

where the factor Zo(L,y) is the partition function Eq. 
(|27| of the free propagation. 



A. Direct integration 

The direct integration of the partition function Eq. (|3|) 
for the random force problem V((j>, x) = f{x)<f>(x) is con- 
veniently done within a Fourier representation. For tech- 
nical convenience we extend the problem to the interval 
[—L,L] and define the anti-symmetric force and displace- 
ment fields f(—x) = —f(x > 0) and </>(—x) = —(f>(x > 0). 
The relevant quantities in Fourier space are the sine- 
transforms 



<]n 



dx g(x) sm(k m x) 



(35) 



with k m = ttitt/L] the back transformation reads 

1 OO 

g(x) = y^ g m {x) sin{k m x). (36) 

m— 1 

The Hamiltonian (on the interval [0, L]) assumes the form 
(we make use of the solution (f> qm = — fm/ck^ of Eq. (|3"Tj) 
in Fourier representation) 



HI 



H.,. 



sr 



J n 



^ 4cLfc 2 



(37) 



y ^ (-i)' 



L *-^ k, r 

m=l r 



' In 



and the partition function reads 

Z(L,y;f) = exp(-^) (38) 

The disorder average in the partition function Eqs. (|3"4")l 
or (|38[) has to be taken over the distribution function for 
the random force /, cf. Eq. (|18|) . or in Fourier space, 



P(fm) = 



1 



V 'AttuL 



exp(-/„/4wi). 



(39) 



Taking the result (j3"5)) to the n-th power and integrat- 
ing over the force distribution Eq. (l39|) , we obtain the 
intermediate result 



Z {L,y) =cxp(^- 

n[ 



xlll 1 

m=l 



2L 

-1/2 



(40) 



exp 



OO 



j3n cy 2 2s 



^ 2L (n 2 m 2 



with 



(3nuL 2 



(41) 



Using the product and partial fraction expansion of cir- 
cular functions 23 



sin • 






OO 

n[ 



i - 



= i+E 



2.s 



tan y/s 
we obtain the final result 
Z n (L,y) = Z(s-e) 

\siii\/sJ 



7U — 1 



(42) 
(43) 

(44) 



1/2 



exp 



S\/S 



tan \/s 



(45) 



with the dimensionless displacement parameter 

= £_JL 

6 2u L 3 " 

With our normalization, the partition sum does not de- 
pend on temperature any more (note that in the calcula- 
tion of the free-energy distribution function, the variable 
s will be integrated over, cf. Eq. (f6"U|)'). The result Eq. 
(l4"4l is well defined provided that < s < it 2 ; the singu- 
larity at s = 7r 2 will determine the shape of the left tail 
in the free-energy distribution function, see (|67| below. 

As a simple application, we can use the partition 
function Eq. (|3"8")) to find the free energy (F) — 



—ThiZ(L,0;f m ) of the string starting and ending in 
= 0. Taking the disorder average over the term 
Z~2 m Pfm/^cLk^ in the logarithm of the partition func- 
tion Eq. (|38p . we obtain the result 



(F) 



■CXI- 



uL 2 



2cn 2 ^— ' m 2 

m— 1 



Uc 

12 



& 



(46) 



where we have used the Riemann Zeta function £(2) = 
7r 2 /6 and the definitions Eqs. (ITT1) and (H"2"|) . Alterna- 
tively, we can use the Eqs. (|16|) and (j44|) and calcu- 
late <F)(L,0) = -(ul 2 /c)a s ^(s;0)| s=0 . With Z(s;0) « 
1 + s/12 we then easily recover the above result. Note 
that the result Eq. (|46|) measures the free energy F with 



respect to the entropic contribution F^ s = T In y/2irLT/c 
of the free string. 

The result for the free boundary condition d x (f>\ x= L = 
is obtained by using an alternative expansion: first, we 
symmetrically extend the system from the interval [0, L] 
to the interval [0, 2L] with the definitions <fi(L + x) = 
4>{L — x), f(L + x) = /(£ — x). Second, we expand 
the analysis to the interval [-2L, 2L] using the same 
anti-symmetric extension as above. As a result, we 
can expand the displacement and force fields into modes 
sm(q m L) with q m = (2m — l)ir/2L, m — 1, . . . , oo and 
hence zero slope at x = L. Following the same steps as 
above, we arrive at Eq. f4"0|) with y = and the product 
corresponding to the expansion of the cosine 2 ^, 



cos yfs = JJ 1 - 

m— 1 



4.s 



7r 2 (2m-l) 2 



(47) 



The final result for the partition function with free 
boundary conditions then reads 



Z(s) = 



1 



a/cOS- 



(48) 



where the regime of applicability is restricted to the do- 
main < s < 7r 2 /4; again, the singularity at s = 7r 2 /4 
determines the shape of the left tail in the free-energy 
distribution function, cf. (|7ip . The alternative procedure 
of realizing the free boundary condition via integration 
over the end-point coordinate y, cf. Eq. ff25]) . provides 
the identical result, although via a much more tedious 
calculation of determinants. 



B. Replica approach 



The replica Hamiltonian Eq. (|20|) for the random force 
problem Eq. © reads 



i-L , n 



1 r t 

HnlM] = 2 / dx[cY,[dMx)] 2 (49) 

n 

-/3w X </> a (x)</)6(a:)| 

a, 6=1 
= -- / fife ^ 4'a{x)U a b4'b{x) (50) 



with the matrix 



C4& = c <5 a6 <9 2 + /3m. 



(51) 



Accounting for the random shift Vq(x), cf. Eq. <JSj> , adds 
an additional term — n 2 /3[/ L/2 to the Hamiltonian (|49|) . 
The matrix U a b can be easily diagonalized and we find 
one (n — l)-fold degenerate eigenvalue Ai = cd 2 per- 
tinent to the free string with the (n — 1) orthonormal 
eigenvectors £f obeying the constraint E"=i £? = 0' 
i = 1, . . . ,?i — 1. The n-th eigenvalue A2 = c<9 2 + (3nu 
is non-degenerate and appertains to an inverted har- 
monic potential problem; the associated eigenvector is 
££ = \j\fn, a = l,...,n. The coefficients £? of the 
(n x n) transformation matrix £? satisfy the conditions 
Ea=i £?£? = <% (completeness) and E™ =1 CCf = <U 
(orthonormality) . In terms of the new fields (pi(x) and 
boundary conditions qi = <Pi(L), 



wo*) = x^-o*)' 



?« = X ^ ay °' 



(52) 



the wave function (or propagator) Eq. (fT9l) takes the form 



*«%};£) =n 



i=l 



D[^(s)] 



2% n (x)] (53) 



x exp 



x exp 



/3c 
/3c 

^"7o 



L n-1 



dx X^ViW] 5 



da: [(^(^(x)) 2 - ip 2 n (x)/A 2 



where we have introduced the length parameter A, cf. Eq. 

(ED, 



f3nu 



__ 

s 



(54) 



The (n — 1) free propagators are given by ^Q^q^L), 
i = 1, . . . , n — 1, cf. Eq. (|2T|) . The propagator \J/ ih (g„; L) 
for the inverted harmonic potential problem is obtained 
by solving the imaginary-time Schrodinger equation 



d x ^ ih (q; x) 



-d 2 + ^q 2 
/3c q A 2 q 



y ih (q;x) (55) 



with the initial condition ^>ih(q', x = 0) = S(q). With the 
Gaussian Ansatz <& ^(q; x) — %(x) exp[— a(x) q 2 /2] and 
proper accounting of the initial condition, we find the 
solution (cf. Ref. [241 ) 



1/2 



/ \/ S \ -V* / 

*ih(q n ;L)=l-rt- r ) exp 
vsin-i/s/ V 



/3c y^s 2 
2L tan ^/i 



(56) 



Inserting the free (Eq. (|2T|) ) and harmonic (Eq. (J56|) ) 
factors into the full propagator Eq. (|53"|) and transforming 



back to original variables, E™=i 9? = E"=i ^f — 
E"=i 2/a - (l/«)(E"=i y°) 2 , we obtain the result 



'17, 



n 

^({y a };L)=\jl^o(ya;L) 



a=l 



gih(EfcWV^£) 



(57) 
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Choosing the appropriate boundary conditions y a = y, 
a = 1, ...,n, we obtain the replica partition function 
Eqs. (fl"9]) and (|2"1~T) identical to the previous result Eq. 
dMl), Z r {s; e) = Z(s; e), < s < it 2 . 

The result for the (replica) partition function has been 
derived for positive integer n. Since Z(s; e) depends on 
n only via the parameter s, the expression Eq. (144)) can 
be analytically continued to the complex half-plane re- 
stricted by the condition Re[s] < ir 2 . 

For the free boundary condition, we obtain the replica 
partition function via integration of Eq. (|57[) over all end- 
points {y a }, cf. Eq. (|26|) : the integration is conveniently 
done in the variables q a and we make use of the nor- 
malization Eq. (j2"9")l to obtain the result identical to Eq. 



|, Z r (s) = 1/y/oos^s, < s < 7r 2 /4. Furthermore, 
the analytic continuation to real negative values of the 
parameter n provides the expression 



Z{s) 



cosh 



(*<0); 



(58) 



alternatively, this result is obtained via the solution of the 
Schrodinger equation Eq. (j55|) for negative n involving a 
summation over the discrete spectrum of the parabolic 
potential, see Appendix |A"1 



C. Distribution function: fixed boundary condition 

We now turn to the calculation of the free-energy dis- 
tribution function VL,y{F) from the partition function 
Z n (L,y). Following the procedure described in Sec. HT1 
specifically Eqs. (jT4"|) and ([To]) , the Laplace transform and 
its inverse assume the form 



Z(s;e) 



+ 00 



dfPe(f) exp(-s/), 



(59) 



i /-Moo 

Pe(f) = 7T- / dsZ(a;e)exp(*/), (60) 



where 



f(F,L) 



F 



FfW 

t T \ c y 



Ff(L) 



-L 2 



(62) 



are the rescaled free energy of the system and the rescaled 
displacement parameter; the original free-energy distri- 
bution function VL,y(F) then derives from the rescaled 
expression p e (f) through the relation 



FlAF) 



1 



Ff(L) 



p t{y ,L)(f(F,L)). 



(63) 



Note, that the parameter £ drops out in the combina- 
tion U c /L 2 , as has to be the case for the random force 
model where the disorder is characterized by only one 



parameter, its strength u. Or in other words, using the 
results below for a random force (rather then a random 
potential) problem, these are valid for all length scales. 
For the relaxational free-energy distribution function of 
the system with fixed boundary condition, we obtain the 
expression 



1 f +lco / 

Pe(f) = 7T- / ds(- 



x exp 



sinys 
tani/s 



1/2 



+ fs 



(64) 



which simplifies drastically for the special case of fixed 
boundary conditions with 0(0) = <j>(L) = 0, 



Po(f) 



1 
2?ri 



ds 



1/2 



(-^ P ) 1/ ex P (/ S ) 
V sin Jsj 



(65) 



The above result already expresses an important prop- 
erty of the distribution function po{f)- for / > the ex- 
pression under the integral is analytic and quickly goes 
to zero at s — > — oo, hence the contour of integration in 
the complex plane can be safely shifted to — oo. This 
implies that the function po(f) must be equal to zero for 
/ > and the relaxational free energy of the directed 
polymer with zero boundary conditions is bounded from 
above, F < 0. This constraint then is easily understood, 
as the presence of a random force can only reduce the 
relaxational free energy of the directed polymer. 




FIG. 2: Relaxational free-energy distribution function p e (f) 
for the randomly forced directed polymer for several val- 
ues of the dimensionless displacement parameter e = 
{y/V) 2 {Lc/Lf/2: e = 0, 0.001, 0.005, 0.01. 



The evaluation of the inverse Laplace transform Eq. 
(|M]) is discussed in Appendix [B] and provides the free- 
energy distribution function p e (/) as illustrated in Fig. [2] 
(note that all temperature dependence has vanished since 
we measure our free energy with respect to the entropic 
contribution / fs = F^/F } = (cT/uL 2 ) In y/2nLT/c of 
the free string). The various asymptotic forms of the 



distribution function p e (f) in the limits / — \ ±00 and 
/ — > —0 are derived via a saddle-point integration and 
read, 



p e (f -> -00) 
p e (f ^+00) 

Pe=0(f "► -0) 



exp 
exp 
exp 



V|/| 



27e 2 W J 

V 32J7|) ; 



(66) 
(67) 
(68) 



note that the shape of the left tail is determined by the 



singularity of Z(s; e) at s 



cf. Eq. flU}. The above 



results agree with those obtained before in Ref. [2C 



D. Distribution function: free boundary conditions 




FIG. 3: Relaxational free-energy distribution function p(f) 
of the randomly forced directed polymer with free boundary 
conditions. 



The result (|48|) provides us with all the moments of the 
relaxational free-energy distribution function, of which 
the first one, the average free energy, is given by 

JF)=-U C (L/L c f [d 8 Z(s)} 8=0 = -^(^) 2 - (69) 

In order to obtain the full distribution function, we per- 
form the inverse Laplace transform (Vl(F) = p(f — 
F/F f (L))/F f (L)) 



P(f) = 7^-1 

2m J_ 



1 



+ 200 

ds — ^=^= 

ioo V C0S V^ 



exp(fs). 



(70) 



Given the scaling form in / = F/Ff(L), the result is 
valid at all scales. Again, for / > the integrand is ana- 
lytic and rapidly approaches zero as s — > — 00 and hence 
the function p(f) must vanish identically for / > 0. The 
functional form for / < is found as before, see Ap- 
pendix |5] The relaxational free-energy distribution func- 
tion V(f) assumes a universal form with no parameters; 
it vanishes identically for / > and its overall form is 
shown in Fig. [3j 

Note that the free energy of the 'trivial' configuration 
4>{x) = is equal to zero and any deviation due to the 
action of the random force can only reduce the energy, 
providing a simple explanation for the cutoff at positive 
energies. The asymptotic behavior in the limits / — » —00 
and / — > —0 can be easily estimated by a saddle-point 
calculation, 



P(f 



-00) 
-0) 



exp 
exp 



(-Tl'l)- 

V 32j/f)- 



(71) 
(72) 



have to account for the random shift Vo(x), cf. Eq. ©. 
Here, we concentrate on the situation with free bound- 
ary conditions. The multiplication of Eq. (|48j) with the 
Gaussian exp((3 2 n 2 UoL/2) — exp[s 2 (L c /L) 3 /2] and sub- 
sequent Laplace transform of 



Z r (s) 



-\/cos- 



1 s^l 

-e 2 ^ 



(73) 



generates the (rescaled) free-energy distribution function 
p t (f) shown in Fig. |4j 




1 / 2 



FIG. 4: Free-energy distribution function p (/) of the ran- 
domly forced directed polymer with free boundary conditions 
including the random shift Vo(x). For L 2> L c , the relax- 
ational part of the free energy dominates the distribution; in 
the regime L < L c , where the shifted random force model pro- 
vides an approximation to the random polymer problem, the 
free-energy distribution function is dominated by the Gaus- 
sian part originating from the shift Vq(x). 



E. Shifted random force model 

In order to find the distribution function for the to- 
tal free energy (rather than its relaxational part), we 



The (total) free-energy distribution function V 1 derives 
from a convolution of the distribution function V^ '(F) of 
the relaxational free energy, cf. Eq. ([70)1 and Fig. [3J and 
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the factor V v originating from the random shift Vo(x), 

/oo 
dF'V f {F')V v (F~F'); (74) 

-oo 

the contribution V from the random shift assumes the 
simple Gaussian form 



v v (f) 



1 



'2nF v {L) 
with the scaling parameter 



eM-{F/F v {L)f/2] (75) 



F V (L) = U cy /L/L c . 



(76) 



For L ^S> L c , the result coincides with that for the re- 
laxational energy, cf. Eq. (|T0[) and Fig. [3l and scales as 
F/Ff(L). In the limit of short lengths L < L c , where the 
model can be used as an approximation of the random 
potential problem, the distribution function is dominated 
by the Gaussian due to the random shift Vq(x). More 
specifically, we find the width of the distribution's body 
through the calculation of the second cumulant: expand- 
ing Z r (s) for small values of s and using Eq. (|16p . we 
obtain 



Z r (s) 

JF) 
F^-F 2 



1 



s s 

+ l + ~2 

~4 
L 
L c 



J_ 
48 



<tJ 



U 2 ^- 



1 ( L 



(77) 
(78) 
(79) 



The leading term ex U c /L 2 — u/c in first moment Eq. 
(I78p derives from the random force part of the disorder. 
This is different from the second cumulant in Eq. (1791) . 
where the first term ex U 2 /L c = Uq derives from the 
random shift Vq(x) and dominates over the contribution 
from the random force (second term ex (u/c) 2 ) at short 
lengths L < L c , hence the width of V 1 is given by Fy 
(see below for a discussion of the corrections ex (L/L c ) 3 ). 
Besides the first two moments/cumulants, we can eas- 
ily determine the scaling of higher moments. Starting 
from the convolution Eq. (|74l) , we note the different scal- 
ing of the arguments in the two distribution functions, 
ex F/L 2 for the relaxational part V* and ex F/y/L for 
the random shift part V V (F). Hence, for small dis- 
tances L, the function V* (F') peaks narrowly near zero, 
while V v (F — F') retains a broader shape; expanding 
V v (F ~ F') around F and integrating over F', we obtain 
the following expansion for the total distribution func- 
tion, 



-pr/ 



V\F) « V V (F)-V V (F)F 



(80) 



-■p v (F)F 2 ' 



where V (F) is the derivative of V with respect to 

/ 

the argument F and (...) denotes averaging over the 



£ 

random force part V* . Using the scaling F ex uL 2 and 

F 2 ex u 2 L A for the moments of the relaxational free- 
energy, we can calculate the dependence of the moments 
F k on the length L. Thereby, we exploit the fact that the 
leading term V v (F) in the expansion (180)) is symmetric in 
F, cf. Eq. (|75p . and hence determines the even moments, 
while the next term is anti-symmetric and generates the 
odd moments; finally, the third term provides the cor- 
rections to the even moments. The combination of the 
scaling of V v (F) (deriving from the random shift Vo(x)) 

and of the first two moments F and F 2 (deriving from 
the random force f(x)) then generates the following non- 
trivial scaling of the moments with different powers in L 
for the even and odd moments, 



p2k 



F 2k ex L k +0(L k+3 ), 

k+2 



ex L K 



(81) 
(82) 



where 0(L k+3 ) denotes a correction term ex L k+3 . In 
particular, F ex uL 2 , f 2 ex UqL + 0(u 2 L a ), F 3 ex uUqL 3 , 
F* ex U§L 2 + O{u 2 U L 5 ), F^ ex uU§L 4 . 

Finally, we can estimate the tails of V t from the convo- 
lution Eq. (1741 using the asymptotic behavior of V^ and 
V v and find a left tail V\F < -F tai i) ex exp(F/F f ) 
and a Gaussian tail on the right, V l (F > Fmn) ex 
exp[-{F/F v ) 2 /2}, where F taii = F v [l + (L c /L) 3 / 2 ]. On 
short scales L < L c , we have i^taii ~ Fy(L c /L) 3 ^ 2 > Fy 
and the random force behavior appears only quite beyond 
the body. 



F. Joint distribution function 

We add a note on the joint free-energy distribution 
function V L , y (F,F'), where F = (F+ + F_)/2 and F' = 
(F + — F-)/2 denote the mean free energy and the free- 
energy difference for two polymer trajectories starting at 
the origin 4> = at x = and ending in the symmetric 
points cf) — ±y at x = L, F± = F(L, ±y; V). Opposite to 
the (^-correlated potential, cf. Ref. [I|| the present case of 
the random force model is less reveiling and we keep the 
discussion short. 

Starting with the original (random) Hamiltonian 
H[<f)(x);V] with the random force potential Eq. ([S]), we 
account for the boundary condition <f>(L) = y through the 
shift <j)(x) — > (y/L)x + <fr(x) (the T = solution for the 
string ending in </>(£) = y derives from the solution end- 
ing in 4>(L) = by adding the shift (y/L)x) and obtain 
the Hamiltonian described by Eqs. (|33)l and (|37|) . The 
relaxational free energy of the system with the boundary 
condition <j>(L) = y separates into the terms 



F[L,yJ} = 



cy 



y 



2L L 



dxxf(x)+F[L,0;f]. (83) 



The first term is the trivial part of the elastic energy, the 
second is a random constant, and finally, the third is the 
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(random) relaxational free energy of the polymer with 
zero boundary conditions; its randomness is correlated 
with the randomness in the second term. Then, for the 
free energies F and F' introduced above, we find that 



F'[L,y;f] 



F[L,y,f\ 



II 



cy 
2L 



dxxf(x), 



F[L,0;f], 



(84) 
(85) 



and hence F' and F carry the information on the second 
and third terms in Eq. (l83l) . respectively. Although the 
joint distribution function for the random and correlated 
quantities F and F' must be non-trivial, we can conclude 
that the separate statistics of F and F' must be simple: 
according to Eq. (|84|) . the distribution for F' is Gaussian 
with zero mean and width (F 1 ) 2 = y 2 uL/3, while the 
distribution for F must coincide with that for the free 
energy with zero boundary conditions Vl, v =o(F), cf. Eq. 
(|65|). shifted by the trivial elastic term cy 2 /2L. Also note, 
that a change in the final coordinate y modifies the poly- 
mer's trajectory over the entire length L and hence the 
joint distribution function is not expected to factorize, 
in contrast to the results found for the short-range cor- 
related random polymer problem^. The detailed replica 
calculation, which represents a straightforward extension 
of the above analysis, produces results in full agreement 
with these simple arguments. 



The propagator for the harmonic-correlator approxima- 
tion then assumes the form (cf. ([57)0 



n 

V({y a };L)=[l[* h (y a -L) 



o=l 






(87) 



where ^h derives from \&ih by the substitution A — > i\. 
For simplicity, we only consider the model with free 
boundary conditions and find the shifted random force 
result Eq. ([73]) replaced by the expression (s = L 2 /X 2 = 
n(U c /T) 2 (L/L c ) 2 ; see also Appendix © 



Z r {s) 



1 



-,(«-!) 



ycosh y/s 



exp[s 2 (L c /L) 3 /2] 



[cOSh y/s\ 2 [cOSh y/s\ 



2U C L* e —J]S 



Although the inverse Laplace transform can be per- 
formed, the resulting (total) free-energy probability dis- 
tribution p(f) develops a negative right tail at zero and 
low temperatures, see Fig. [5] at large temperatures T » 
U c the right tail exhibits pronounced oscillations. These 
unphysical results are due to the departure of the ap- 
proximate harmonic interaction U p (<p) from the true in- 
teraction U ((/)), becoming relevant at large scales L > L c , 
<f> > £, F > U c , and the large fluctuations of the string 
at high temperatures T 3> U c . Note that the inverse 
Laplace transform cannot be performed at all in case the 
random shift Vq(x) is ignored. 



IV. HARMONIC-CORRELATOR 
APPROXIMATION 



We consider the random directed polymer described 
by the Hamiltonian Eq. (jTJ) and approximate the inter- 
action —f3 2 U in the replica Hamiltonian Eq. (|20[) by the 
harmonic expression Eq. ([9]) to arrive at, 



H n [{4>a}} = / dx\-Y^[d x M^)\' 

J I Z a =l 



(86) 



0U 



-/3U L 



a,b=l ' 

1 f L - ~ n 2 

~x j dx ^2 ^(x) U ab c/) b (x) - —f3U L 

"'° a,b=l 



with the matrix U a b = (cd 2 — (3nu)6 a b + /3u (note that the 
parabolic approximation of the correlator should be im- 
plemented after the integration over the disorder poten- 
tial). Diagonalization produces the (n — l)-fold degener- 
ate eigenvalue Ai = cd 2 — j3nu of the harmonic oscillator 
problem with the (n— 1) orthonormal eigenvectors £° con- 
strained by the condition X)"=i ^f 1 ^ 0, i ^ 1,...,/t, — 1, 
and one non-degenerate eigenvalue A2 = cd 2 of the free 
problem with the eigenvector £° = 1/y/n, a = 1, . . . ,n. 




f 5 



FIG. 5: Free-energy distribution function p(f) for the directed 
polymer with free boundary conditions using the harmonic- 
correlator approximation. Solid curves refer to T — 0, while 
the dashed curve attains to T = 2U C - 



The breakdown of the harmonic-correlator approxima- 
tion is conveniently observed in the second moment: ex- 
panding Z r (s) for small values of s we find to second 
order 



Z r {s) ~ 1 



o2 r 



1 



T L 

W C T C 



1 



(89) 
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and using Eq. (TT6l) , we find the average free energy F = 
-U C (L/L C ) 2 '/A (cf. Eq. ([75]) ) and the second cumulant 
reads 



f^-f 2 = u 2 — 



T L 

2U~ C L, 






1 (L\- 
12 VZ 



(90) 



Comparing with the result Eq. (|79|) for the shifted ran- 
dom force, we note the additional dependence on tem- 
perature and the sign-change in the correction term 
oc (Lj L c ) 3 . This decrease in width is in accord with the 
behavior of the free energy fluctuations in the random di- 
rected polymer problem ([2]), which scales as SF oc L 2 ^ 1 
at large distances; with a wandering exponent 6 £ = 2/3, 
we have SF oc L 1 / 3 . The negative correction ex (L/L c ) 3 
to the linear growth in 5F 2 observed in Eq. (1UU1) then is 
consistent with the sublinear growth 8F 2 ex L 2 ^ 3 of the 
exact solution. 

As before, we can analyze the higher moments of the 
distribution function and compare the results for the 
harmonic-correlator approximation at T = with those 
obtained for the shifted random force model. Making 
use of the product form of Z r and expanding each fac- 
tor in s, we find identical leading terms for all even 
and odd moments (corresponding to equal expressions 
for the first two terms in Eq. ([80]) ). while the correc- 
tions to the even moments (described by the third term 
in Eq. ([80]) ) are different. Furthermore, we note that 
the even moments F 2 * oc U 2k {L/L c ) k [l - 0{L 3 /L 3 c )] oc 
U§L k are large with a small negative correction, while 
the odd moments F 2k + l oc -U 2k+1 (L/L c ) k+2 [l + 
0{L 3 /L 3 )] oc uUQL k+2 are small, their ratios being 
(F 2k+1 ) 2 /F 2 ( 2k + 1 ') oc (L/L c ) 3 . To leading order, the 
free-energy distribution function for the random poten- 
tial model at small scales then is a trivial Gaussian gen- 
erated by Vo(x), with a small negative shift and a small 
reduction in width due to the random force term in the 
potential. 

At T - 0, the second cumulant turns negative for 
L > 2^/3/2 L c and the result Eq. (|90|) makes no longer 
any sense, hence the harmonic approximation to a ran- 
dom potential problem cannot be used on scales larger 
then £ (along the transverse direction) or L c (along the 
longitudinal direction); at finite temperatures the regime 
of validity is further reduced. 

Although the results for the shifted random force re- 
main valid at any length L, we emphasize that the har- 
monic correlator provides a better approximation for the 
behavior of the short-range correlated random polymer: 
Both results agree in lowest order, providing the same 
first moment F due to the random force f(x) and the 

same leading term in the second cumulant F 2 — F 
generated by the random shift Vq(x). The correction 
oc (LjL c ) 3 in the second cumulant is due to the ran- 
dom force f{x) and contributes with the opposite sign 
in the shifted random force model as compared to the 
harmonic correlator approximation. While the shifted 
random force result Eq. (j?9")) is correct (at all scales) 



when dealing with a true random force model, the correc- 
tion oc (L/L c ) 3 carries the wrong sign when used as an 
approximation to the random potential model and it is 
the result of the harmonic correlator approximation Eq. 
(|90| which should be trusted. Indeed, the harmonic cor- 
relator preserves the translation invariance of the prob- 
lem, whereas some quadratic terms are dropped from the 
shifted random force model. Expanding the potential to 
to second order, 

V(cf>, x) = V (x) + f(x)</>{x) - \g{x)<p(x), (91) 



we identify the terms in V(</>, x)V(cj>',x') with the har- 
monic expansion Eq. © to obtain the correlators Eqs. 
(J7J and ((8j, Vo(x)g(x') — uS(x—x'), and vanishing mixed 
terms Vo(x)f(x') = f(x)g(x') = 0. A scaling estimate of 
the second moment F 2 using Eq. (I^Tjl then provides a 
leading term oc L from Vq and subleading terms oc L 4 
from / and from g. The contribution from the random 
force provides the positive contribution U C (L/L C ) 4 / 12 in 
the cumulant Eq. (|79l) . while the mixed terms Vo(x)g(x') 
contribute with a negative weight, generating the nega- 
tive correction oc L A in Eq. (|90|) . Note that higher order 
terms do not change this result but contribute to the next 
order term oc L 7 . 

Given that the harmonic correlator provides the better 
approximation to the random polymer problem at short 
scales, one may wonder why we end up with unphysical 
results (negative distribution function, negative second 
moment) at larger scales. Also, different types of cor- 
relators, e.g., power-law type, have been studied in the 
past, cf. Refs. [0 and [2l|, and one would like to know, 
what properties of a disorder correlator guarantee con- 
sistent results; this question is addressed in the following 
section. 



A. Correlators with non-positive spectrum 

It is important to identify problematic correlators right 
from the beginning; indeed, the proper definition of the 
disorder potential is subject to important constraints 25 
regarding its shape U{4>) and failure to respect these 
constraints may lead to unphysical results. Consider a 
random potential V(</>) and its Fourier representation 
V(q) = J d(f>V(cj)) exp(— iq<f>); then the Gaussian distribu- 
tion function of the random function V(q) has the form 



-dm mp ( f d i lis, n 

nv( q )} = Poe^-j-^^y, 



(92) 



the width G(q) has to be positive and relates to the cor- 
relation function U(<p) via 



^) = |^G(9)exp(i#), 



(93) 
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Expanding both sides in powers of < 



u(o) + J2 



t/( 2fc )(0) 



(94) 



fc=i 



fe=i 



(~l) fc 

(2fc)! 



/|«(^ 



2fc\ /2fc 



and comparing coefBcients, we find that the 2fc-th 
derivative of U(4>) in the origin relates to the integral 
J dqG(q) q 2k , which is a positive quantity. Hence, we 
have to be careful in our choice of the correlator U(cj>). 
E.g., truncating the expansion of U((f>) beyond some fc*, 
such that U^ 2k \0) = for k > k* , we impose the condi- 
tion 



f dqG(q)q 2k = for k > k* 



(95) 



which cannot be satisfied for a positively defined G(q). 

Obviously then, choosing a parabolic correlator U p ((j>) 
as in Eq. ([9]) is in severe conflict with the constraint Eq. 
([95]) . The averaging over the disorder potential V{<j>, x) 
is undefined for those modes (in Fourier space) where 
G(q) is negative. Hence, going over from the disordered 
directed polymer (a statistical mechanics problem) to 
the quantum boson problem is an ill-defined step and 
the results cannot be trusted any longer. On the other 
hand, performing the integration over the random poten- 
tial V(<p,x) with a well defined, i.e., spectral, correlator 
U{4>) and expanding the resulting interaction —/3 2 U((j)) 
between bosons is perfectly admissible and leads to an 
identical result; in this case, however, we know that 
the quantum boson problem does not describe the ran- 
dom polymer problem on scales where the approximate 
quadratic correlator deviates strongly from the original 
correlator. Nevertheless, in the end we have to appre- 
ciate, that the harmonic-correlator approximation ([9]), 
although breaking down at lengths beyond L Cl does pro- 
duce more accurate approximate results for the short- 
range correlated random potential problem (|5J| than the 
shifted random force model ([5]), although the latter re- 
mains formally valid at all length scales L. The (shifted) 
random force model then should be used whenever the 
disorder landscape is given by a force field as defined by 
Eqs. (0, 0, and ©, but not as an approximation to a 
random potential problem. 



(y) 2 are easily calculated within replica theory 21 . Defin- 
ing 



(y a yb) 



n „ 

n m* 

0=1*' 



y a yb^({y c };L) 



we obtain the two types of averages 



(y 2 ) = lim(y a y b )\ a=b , 

n—tO 



(96) 



(97) 



(y) 2 = lim(y a y b )\ a7 t b . 

n— 7-0 



The Hamiltonians for the shifted random force model 
and the harmonic-correlator approximation differ only by 
the term (/3un/2) X)"=i ^a> which vanishes in the limit 
n — > 0, hence both schemes produce identical results 
for the displacement correlators in Eq. (|97|) . We then 
concentrate on the random force case and calculate the 
expression 



(VaVb) = C 



II / dy c y a y b exp[--^2K cd y c y d ] (98) 

c=l"' J c,d 



with K r 



Ad, 



cd 



B and 



A- £ 

A ~ L> 



nL 



C 



/ /3c y /2 ( Vs \ 
\2ttL) \sin,/s) 



tan y/s 

1/2 



- 1 



In the calculation of (y a yb)\a^b, we combine all diagonal 
terms into a sum (D/2) J2 C Vc with D = A + B, leaving 
the non-diagonal in the form (5/2) ^2 c ^ d y c yd\ the non- 
diagonal average then follows from the derivative 



(y a yb)\ a ^b 



d 



n „ 

Ujdy c 



*({y c };L)\ D , 



n(n — 1) dB 
while the diagonal average is given by the derivative 



(yl) 



2_d_ 

ndA 



n „ 



n{yc};L)\ 



(99) 



The final results assume the form 



B. Displacement correlator 

Another quantity of interest in the random poly- 
mer probem is the displacement correlator ((j> 2 )(L) = 
{[(j){L) - <p(0)} 2 ) , with (...) and (. . . ) denoting ther- 
mal and disorder averages, respectively. Choosing free 
boundary conditions with 0(0) = and an arbitrary po- 
sition <p(L) = y for the end-point, the averages (y 2 ) and 



W) -H»fi-^ -fL± + £&* 



(y) 2 = 



n-*)\A A 2 J 
— lim 



o(D-B) 



U c L t 



lim — 
n-X) A 2 



4(r).« 



3 \L B 



the relation (y 2 ) — (y) 2 = (y 2 )\v=o = TL/c (here, 
{y 2 )\v=o denotes the thermal average in the absence of 
any disorder, V = 0), is a constraint holding true for any 
disorder potential uncorrelated in x, cf. Refs. [26J andl27l. 
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V. SUMMARY AND CONCLUSIONS 

The (shifted) randomly forced polymer model and 
the random disordered polymer described through 
a harmonic-correlator approximation define quadratic 
problems and hence admit exact solutions. For the ran- 
dom force models, different approaches can be taken, ei- 
ther a direct integration of the path-integrals within a 
Fourier representation or using the (real space) replica 
technique; in retrospect, the preferred method is a matter 
of taste. We have determined the free-energy distribution 
functions Vl, v (F) and Pl{F) for fixed and free boundary 
conditions. This calculation ne cessitates th e determina- 
tion of all powers Z r (n;L,y) = Z n [L,y;V] (rather than 
the usual n — > limit) and subsequent inverse Laplace 
transformation of the analytically continued replica par- 
tition function Z r (n e C; L,y). The displacement corre- 
lators (y 2 ) and (y) 2 have been found as well. The sim- 
plicity of the quadratic models allows to carry through 
the entire program and thus serves to study not only the 
physical properties of the problem but its methodological 
aspects as well. 

Regarding the shape of the distribution functions 
for the random force model, a number of interesting 
features has been obtained: for the free boundary, 
the probability to find a positive free energy F van- 
ishes exactly, with an essential singularity appearing 
in Vl(F) oc exp(— uL 2 /32c\F\) as F approaches zero 
from the left, cf. Eq. [72] For fixed boundary condi- 
tions, a similar result has been found for Vl, v =o{F) 1 
see also Ref. [20|. Furthermore, the left and right tails 
provide a consistent scaling F <x L 2 and y oc L 3 / 2 , 
V L , y {F -> -co) ex exp(-7r 2 c|.F|/wL 2 ) and V L . y (F -*■ 
oo)'oc exp(-(16/27)|F| 3 /ucy 4 ), cf. Eqs. §6§ and (joTjl . 

When interested in the short distance behavior of the 
random directed polymer Eq. (J2J, two types of approx- 
imations offer a drastic simplification of the problem: 
these are the expansion of the random potential V((f>, x) 
according to Eq. ([5]) (generating the shifted random force 
problem) or the expansion of the correlator Eq. ([9]) (lead- 
ing to the harmonic-correlator approximation). While 
both approximations generate the same results for the 
even and odd moments to leading order, the next to lead- 
ing order terms turn out different. In this situation, the 
results of the harmonic-correlator approximation have to 
be trusted, as it consistently accounts for the relevant 
terms preserving the translation invariance of the prob- 
lem. Collecting all results, we find that the free-energy 
distribution function for the random potential model at 
(T = and) small scales is a trivial Gaussian of width 
U C \JL/L C generated by Vo(x), with a small negative shift 
-U C {L/2L C ) 2 and a small reduction ~(U C /2A)(L/L C ) 7 / 2 
in width due to the random force term in the potential. 

Finally, we mention a few useful insights regarding the 
replica technique which derive from our analysis above. 
The replica technique provides a link between two seem- 
ingly unrelated problems, the classical statistical theory 



of disordered polymers and the quantum many-body the- 
ory of attractive bosons. Several stumbling blocks can 
be eliminated by properly appreciating the subtleties in 
this mapping. As is well known, after the mapping from 
polymers to bosons the disorder correlator assumes the 
role of the interaction potential. While many shapes 
for the interaction potential may produce meaningful re- 
sults for the quantum boson problem, only a restricted 
set of them (those describing a correlator with positive 
spectrum) relate to a meaningful random polymer prob- 
lem. Hence, the original choice of physical correlators and 
any modification thereof during the calculation should be 
done with great care; in particular, a simple power-law 
form 1 might not work. E.g., there is nothing wrong in 
studying quantum bosons with a simple harmonic inter- 
action U(<f>) = —Uo + u(f) 2 /2 and the results obtained for 
the quantum propagator are perfectly acceptable for any 
constant shift Uq. However, interpreting the result for 
the propagator in terms of a replica partition function 
and transforming back (via the inverse Laplace transfor- 
mation) to random polymers, the resulting distribution 
function becomes unphysical when setting Uq = 0; drop- 
ping a shift Uq in the potential for the bosons is a triv- 
ial shift in energy, while ignoring the same shift in the 
correlator produces unphysical results for the polymer 
problem after Laplace transformation. 
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Appendix A: Negative replica number 

We determine the replica partition function Z r (s), Eq. 
(|58[) . for the polymer with free boundary conditions via 
direct solution of the Schrodinger equation Eq. (|55|) for 
negative n. We confirm, that the result analytically 
continued from positive n agrees with the one obtained 
for negative n. The wave function ^f(q,x) satisfies the 
Schrodinger equation (cf. Eq. (|55j) ) 



d x V(q;x) 



1 



1 2 f3c 2 



*(q;x), (Al) 



with A 2 = c//3\n\u. We are seaking the solution 



*(g; x) = J2 A k eM-Ekx) * fc (g), (A2) 

satisfying the initial condition ^(q;x = 0) = S(q); the 
energies and corresponding orthonormal eigenfunctions 
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Ek and ^k(q) satisfy the stationary equation, 



Ek^k(q) = -\ [jf q + ^ q 2 } **(?)• (A3) 



The coefficients Ak in Eq. (|A2[) derive from the initial 
condition 



A h 



+ 00 



dq%(q)y(q;x = 0) = %(0). 



The spectrum of the harmonic problem is given by Ek = 
(k + l/2)/A and the corresponding eigenfunctions are 
(see, e.g.. Ref.[H) 

MQ) = (^^ I ) 1/2 exp[-(/3c/2A)g 2 ]iJ fc [^ c 7A (Z ], 

where Hk(x) are the Hermite polynomials Hk(x) — 
(-l) fc cxp(x 2 )9^[exp(-x 2 )]. Substituting A k and ^k 
into Eq. (IA2|) and taking into account that i?2/+i(0) = 0, 
we obtain the wave function 



In order to find the explicit form of the function R(r]), 
we implement the shift I — > I + 1 in the sum (|A8[) and 
obtain, 



^ ( 2 ;)! , ~ 

R w = i+EkU = 1 +E 



(21 + 2)! 



J+i 



1^2 



^ my 

1=1 y ' 1=0 



7'/' 



/ + 1 



-(/!) 2 ' ^(j + i)(/!)2 
l + 4r]R(r))-2S(r)), (A9) 

(A10) 



^^Etttt^'''" 1 



;=n 



(i + i)(i\y 



With i? the derivative of S, R(rj) = d v S(r]), we obtain 
the differential equation d v S(r]) = 1 + 4r]d n S(r]) — 2S{rf) 
and the initial condition 5(0) = determines the solution 
S(rj) = (1 - VI -4t7)/2, from which #(77) = 1/a/I - 4?; 
follows via integration. Substitution into Eq. (|A7|) pro- 



duces the final result Z(s) — 1/ycosh-^/i, in agreement 
withEq. ([58]). 



*(g;s) 



E 



1 



7rA^^2 2 '(2n! 
z=o v ' 



exp(-i?2;x) 



(A4) 



x exp[-(/3c/A) g 2 ] H^y/fh/Xq] H 2l (0). 

With the spectrum Ek(X) and the normalization 
H2i(0) — (— 1) (2I)!/2'Z!, we obtain the replica partition 
function for free boundary conditions (cf. Eq. (J48J1) 



/+00 
dg*(q;a; = L) 



(A5) 



exp(-V5/2) ~ (-1)' 

Z^^r7r ex P(- 2 vsO^ 



v^ 



/=() 



2 2/ Z! 



where s = L 2 /X 2 and 



Cfc 



+ 00 



dx exp(— x 2 /2) Hk(x). 



(A6) 



Using the recurrence relation Hk+\(x) = 2xHk{x) — 
2kHk-i(x), we find that Ck+2 = 2(fc + l)Cfc and with 
Co = V27T, we obtain the coefficients Ck — \/2tt{21)\/1\. 
Substitution into Eq. (|A5j) provides the replica partition 
function in the form 



Z(s) = V2 cxp(-V^/2) R[v(s)} 
with the function R(r]) defined by the series 



(A7) 



1=0 



(I 



na 



(A8) 



and we have introduced the shorthand r](s) 
- exp(-2v^)/4. 



Appendix B: Inverse Laplace transformations 



The inverse Laplace transforms Eqs. (|64|) and (f70[) are 
reduced to the following expressions: Using the transfor- 
mation s — pexp(±iir/2) in Eq. (I64p . we analytically con- 
tinue the expression for the distribution function p e (f) to 
the imaginary axis, 



1 f°° 
Pe (f) = -Re / dp 
it Jo 



x exp 



v / pcxp(z7r/4) 
Vsin[(l + i)V/>72] 



1/2 



(Bl) 



'tan{(l + i)y/p/2] 



ifP 



A change of the integration variable p = 2i 2 provides, 
after some algebra, the final expression 

o5/2 /-oo 

Pe(f) = / d« 3/2 exp[- £W _(i)] (B2) 

cos[7ft)/2 + 2t 2 / + tt/8 - eo; + ft)] 

x v^y ' 

The functions $(£), w±(i), and 7(f) are defined as, 



$(*) = \/[sin(i)cosh(i)] 2 + [cos(i) sinh(i)] 2 , 

w±(*) = t 3 [sinh(2t)±sin(2t)]/$ 2 (i), 

sin(7(t)) = -cos(*)sinh(*) /$(*), 

cos(7(i)) = sin(i)cosh(£)/$(i). 

Similarly, substituting s = 2i 2 exp(±i7r/2) in Eq. ffTO]) . 
one obtains 



4 f 00 
P / = - dt 
* Jo 



cos[C(t)/2 + 2/t 2 



(B3) 



16 



with the functions $?(t) and £(t) defined by 



*(i) = \/[cos(t) cosh(i)] 2 + [sin(i) sinh(t)] 2 , 
sin(C(t)) = sin(i)sinh(t)/*(i), 
cos(C(*)) = cos(t)cosh(t)/*(t). 



The remaining integrals in Eqs. (|B2[) and (IB3|) have to 
be done numerically. 



1 T. Halpin-Healy and Y-C. Zhang, Phys. Rep. 254, 215 
(1995). 

2 S. Lemerle, J. Ferre, C. Chappert, V. Mathet, T. Gia- 
marchi, and P. Le Doussal, Phys. Rev. Lett. 80, 849 (1998). 

3 G. Blatter, M.V. Feigel'man, V.B. Geshkenbein, A.I. 
Larkin, and V.M. Vinokur, Rev. Mod. Phys. 66, 1125 
(1994). 

4 D. Wilkinson and J.F. Willemsen, J. Phys. A 16, 3365 
(1983). 

5 J.M. Burgers, The Nonlinear Diffusion Equation (Reidel, 
Dordrecht, 1974). 

6 D.A. Huse and C.L. Henley, Phys. Rev. Lett. 54, 2708 
(1985); M. Kardar and Y-C. Zhang, Phys. Rev. Lett. 58, 
2087 (1987). 

7 M. Kardar, Nucl. Phys. B 290, 582 (1987). 

8 Y.-C. Zhang, Europhys. Lett. 9, 113 (1989). 

9 M. Prahofer and H. Spohn, Phys. Rev. Lett. 84, 4882 
(2000). 

10 I.V. Kolokolov and S.E. Korshunov, Phys. Rev. B 75, 
140201 (R) (2007); I.V. Kolokolov and S.E. Korshunov, 
Phys. Rev. B 78, 024206 (2008); see also I.V. Kolokolov 
and S.E. Korshunov, Phys. Rev. E 80, 031107 (2009), 
where the tails in the free-energy distribution are calcu- 
lated for the random polymer problem with a finite-width 
correlator. 

11 V. Dotsenko, 'Bethe Ansatz derivation of the Tracy- 
Widom distribution for one-dimensional directed poly- 
mers', arXiv: 1003.4899 and arXiv: 1004.4455 (2010). 

12 P. Calabrese, P. Le Doussal, A. Rosso, 'Free-energy dis- 
tribution of the directed polymer at high temperature', 
larXiv: 1 002.4560 (2010). 

13 T. Sasamoto, H. Spohn, 'Universality of the one- 



dimensional KPZ equation', arXiv: 1002. 1883 (2010) and 
'Exact height distributions for the KPZ equation with nar- 
row wedge initial condition', arXiv: 1002. 1879 (2010). 

14 G. Amir, I. Corwin, J. Quastel, 'Probability Distribution 
of the Free Energy of the Continuum Directed Random 
Polymer in 1 + 1 dimensions'. larXiv: 1 003.0443 (2010). 

15 D.A. Huse, C.L. Henley, and D.S. Fisher, Phys. Rev. Lett. 
55, 2924 (1985). 

16 G. Parisi, J. Phys. France 51, 1595 (1990). 

17 M. Prahofer and H. Spohn, J. Stat. Phys. 108, 1071 (2002). 

18 V.S. Dotsenko, L.B. Ioffe, V.B. Geshkenbein, S.E. Ko- 
rshunov and G. Blatter, Phys. Rev. Lett. 100, 050601 
(2008). 

19 A.I. Larkin, Zh. Eksp. Teor. Fiz. 58, 1466, 1970. [Sov. 
Phys. JETP 31, 784, 1970]. 

20 D.A. Gorokhov and G. Blatter, Phys. Rev. Lett. 82, 2705 
(1999). 

21 G. Parisi, Rend. Fis. Ace. Lincei IX-1, 277 (1990). 

22 D.A. Gorokhov, Dynamics of Disordered Systems, PhD 
Thesis, ETH Zurich, Diss. ETH No. 13070. 

23 Handbook of Mathematical Functions, edited by M. 
Abramowitz and LA. Stegun (Dover, New York, 1965). 

24 J. Zinn- Justin, Quantum Field Theory and Critical Phe- 
nomena (Clarendon, Oxford, 1993). 

25 D.A. Gorokhov and G. Blatter, Phys. Rev. B 59, 32 (1999). 

26 U. Schulz, J. Villain, E. Brezin, and H. Orland, J. Stat. 
Phys. 51, 1 (1988). 

27 S.E. Koshunov, Phys. Rev. B 63, 174514 (2001). 

28 L.D. Landau and E.M. Lifshitz, in Quantum Mechanics, 
Vol. 3 of Course in Theoretical Physics (Pergamon Press, 
London/Paris, 1958). 



